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In this review we discuss the connection between two seemingly disparate topics, macroscopic gravity on 
'astrophysical scales and Hamiltonians that are not Hermitian but PT symmetric on microscopic ones. In 
i-^ ^particular we show that the quantum-mechanical unitarity problem of the fourth-order derivative conformal 
^ gravity theory is resolved by recognizing that the scalar product appropriate to the theory is not the Dirac 
<D norm associated with a Hermitian Hamiltonian but is instead the norm associated with a non-Hermitian 
i Ibut PT-symmetric Hamiltonian. Moreover, the fourth-order theory Hamiltonian is not only not Hermitian, 
it is not even diagonalizable, being of Jordan-block form. With PT symmetry we establish that conformal 
^ 'gravity is consistent at the quantum-mechanical level. In consequence, we can apply the theory to data, 
tS- 'to find that the theory is capable of naturally accounting for the systematics of the rotation curves of a 
^ -large and varied sample of 138 spiral galaxies without any need for dark matter. The success of the fits 
provides evidence for the relevance of non-diagonalizable but PT-symmetric Hamiltonians to physics. 

in : 
o ■ 

<N ;1 Introduction and Overview 

.£h 'On January 19, 2000 a paper by myself and Aharon Davidson entitled "Fourth order theories without 
^ ighosts" appeared on the arXiv pQ. The abstract of the 2000 paper read: 

"Using the Dirac constraint method we show that the pure fourth-order Pais-Uhlenbeck oscillator model 
" is free of observable negative norm states. Even though such ghosts do appear when the fourth order theory 
is coupled to a second order one, the limit in which the second order action is switched off is found to be a 
highly singular one in which these states move off shell. Given this result, construction of a fully unitary, 
renormalizable, gravitational theory based on a purely fourth order action in 4 dimensions now appears 
feasible." 

Also in the 2000 paper we identified a particular Hamiltonian H(e = 0) = 87a; 4 (2b^b + a^b + b^a) + ui and 
noted of it: "we see that the complete spectrum of eigenstates of H(e = 0)) is the set of all states 
a spectrum whose dimensionality is that of a one rather than a two-dimensional harmonic oscillator, even 
while the complete Fock space has the dimensionality of the two-dimensional oscillator." 

To be specific, what we had found in this paper and its follow-ups [21 [3] was that in the one-particle 
sector the Hamiltonian of the pure fourth-order Pais-Uhlenbeck oscillator model took the generic form 

M=(J J), (1) 
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and that even though the secular equation \M — \I\ = had two solutions (both equal to 1), M only had 
one eigenvector since the eigenvalue equation 
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only had one eigenvector solution, viz. the q = one that obeyed 
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(3) 



The one-particle sector Hamiltonian was thus of a non-diagonalizable, Jordan-block form. And with 
its missing an eigenvector, its eigenvector basis was incomplete (in consequence of which the two energy 
eigenvalues had to be equal since they had to share a common eigenvector), and M could not be diagonalized 
by a similarity transformation. Even though the Hamiltonian was not Hermitian, nonetheless all of its 
eigenvalues (as well as those in the multi-particle sectors of the theory) were real. The Pais-Uhlenbeck 
Hamiltonian thus serves as an example of a non-Hermitian Hamiltonian that has all eigenvalues real, to 
thus show that while Hermiticity of an operator is sufficient to yield real eigenvalues, it is not necessary. 

Aharon Davidson and I had been exploring the Pais-Uhlenbeck oscillator model because it serves as a 
prototype for the fourth-order conformal gravity theory that I had been advancing (see e.g. [HE]) as a 
candidate alternative to standard gravity and string theory. My interest in the conformal gravity theory 
was triggered by the fact that in a conformal invariant theory the cosmological constant would have to be 
zero, with conformal gravity thus serving as a good starting point for addressing the cosmological constant 
problem. Problematically though, when quantized, fourth-order theories such as conformal gravity were 
thought to possess negative norm states and not be unitary. However, through our work Aharon Davidson 
and I were able to show that the theory could solve its ghost problem by having the ghosts not appear 
as on-shell energy eigenstates, to thereby potentially open the road to the construction of a consistent 
quantum gravitational theory in four spacetime dimensions. 

As such, the pure fourth-order Pais-Uhlenbeck oscillator model provides an example of a PT-symmetric 
Hamiltonian at a so-called exceptional point - a point at which the energy eigenspectrum is incomplete 
and the Hamiltonian cannot be brought to a diagonal form by a similarity transformation. However, back 
in 2000 I had no idea that the theory was a PT theory, let alone one at an exceptional point, and only 
learned of it from Carl Bender in 2007 when the two of us started to collaborate on the theory [6J [8]. 
With Carl Bender and his collaborators only having started to develop the concept of non-Hermitian but 
PT-symmetric Hamiltonians in 1998 by showing that such Hamiltonians could have real eigenvalues (see 
e.g. [9j), my 2000 paper with Aharon Davidson thus provides one of the first examples of a theory with a 
non-Hermitian but PT-symmetric Hamiltonian at an exceptional point. 

As first presented by Pais and Uhlenbeck [To], the Pais-Uhlenbeck model is based on an action 
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dt[z 2 - (ujf + w$)zr + uj{uj 2 z 2 } 



(4) 



that describes the dynamics of a single, non-relativistic coordinate z(t). The associated Euler-Lagrange 
equation of motion is of the fourth-order derivative form 
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z(t) = 0, 



(5) 
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and the associated propagator is of the form 



£>pu(w 2 ) 
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(6) 




uj) (u 2 - w|) (u 2 -ul) I (u 2 - uf) (u 2 - Ujf) 



viz. a propagator whose poles are entirely on the real energy axis. 

Two things in the structure of the propagator D-p-^iu: 2 ) immediately stand out: the presence of the 
relative minus sign in the unequal-frequency partial fraction decomposition and the singularity of the 
prefactor in this decomposition in the equal-frequency limit. We thus recognize two distinct sectors of the 
theory, the non-singular case of unequal frequencies (u\ ^ U2) and the singular case of equal frequencies 
(u>i = W2). Prior to my work with Carl Bender the presence of the relative minus sign in the unequal- 
frequency propagator had always been taken to indicate that some of the states in the theory would have 
to have negative norm (since some of the poles would have negative residues), with the theory then not 
being unitary. And the result obtained by myself and Aharon Davidson was that even if there were such 
negative norm states in the unequal-frequency theory, because of the singular nature of the equal-frequency 
limit, in the limit the Hamiltonian would lose some of eigenstates, and none of the energy eigenstates that 
remained would have negative norm. 

My collaboration with Carl Bender then went further by establishing that the unequal-frequency theory 
itself had a realization that was unitary, and that it too could be reinterpreted as a PT rather than a 
Hermitian theory, and that then even in the unequal-frequency case there would be no states of negative 
norm. Thus for both the unequal-frequency and equal-frequency theories the Hamiltonian would not be 
Hermitian. However, since all the poles of the Dpu(u> 2 ) propagator lie on the real axis, all energy eigenvalues 
would be real, with the Hamiltonian then being PT symmetric. At the time Carl Bender and I first studied 
the Pais-Uhlenbeck model we showed that the Hamiltonian was PT symmetric simply by examining the 
behavior of the fields in it under a PT transformation. Subsequently, we established a general theorem 
[TT] that stated that if a Hamiltonian did not have a PT symmetry (as understood in general to mean 
the product of discrete linear and discrete anti-linear operators), some of its eigenvalues would have to be 
complex, with no non-PT-symmetric Hamiltonian then being able to have a set of eigenvalues that was 
entirely real. Since the energy eigenvalues of the Pais-Uhlenbeck theory were all real and the Hamiltonian 
was not Hermitian, the theory thus had to be a PT theory. 

While my own interest had primarily been in the equal-frequency theory since it is a non-relativistic 
prototype of pure fourth-order gravity, the unequal-frequency theory is of interest not only because it is 
the quantum-mechanical limit of the often studied second-order plus fourth-order gravity theory, but also 
because it allows one to monitor what precisely happens in the limit in which the second-order piece is 
switched off. If one ignores metric indices (which as we show below we can), the issues involved can be 
illustrated by consideration of the second-order plus fourth-order scalar field theory with action, equation 
of motion and propagator of the form 




(7) 



(d 2 -V 2 )(d 2 -V 2 + M 2 )#M) = 0, 




(9) 
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where k 2 = —k^ + k 2 . For this propagator a standard Feynman contour integration yields a positive 
frequency contribution of the form 



' j P rp ' - rp rp 



(10) 



where ip n {x) = e lkn ' x , ip m {x) = e lkm ' x , E n = \k n \, E m = (k^ + M 2 ) 1 / 2 . An identification of the propaga- 
tor with the Green's function (Q\T(<j)(x)(j)(x'))\Q) would then require the closure relation for the energy 
eigenstates to be of the form 

— V] \m)(m\ = I, (11) 

n m 

with the presence of negative residues in the propagator translating into the presence of negative norm 
states in the Hilbert space. 

To avoid this outcome we must not identify the propagator with (fl\T((f)(x)<p(x , ))\Q) but must instead 
seek some other identification for it, one in which a closure relation involving negative norm states would 
not then appear. That we might be able to do this at all is due to the fact that the propagator Ds(k 2 ) is a 
c-number quantity that is constructed from a c-number differential equation of the form (d 2 — V 2 )(<9 2 — V 2 + 
M 2 )Ds(x, x') = 5 A (x — x'), and as such it itself makes no reference to any Hilbert space. There is thus some 
freedom in relating the c-number propagator to an underlying q-number field theory, with other possible 
quantum-field-theoretic representations of the propagator being possible. And indeed, it is only after 
second-quantizing the scalar field in a very specific way, viz. the canonical <fi{x) = ^Z{a{k)e lk ' x + a 1 (k)e ' x ) 
where the operator associated with the negative-frequency sector is expressly the Hermitian conjugate of 
the operator a(k) associated with the positive-frequency sector, that we then construct the propagator as 
a matrix element of products of field operators in a vacuum state \Q) that a(k) annihilates on the right 
according to a(k)\Q) = and at(k) annihilates on the left according to (f2|at(A;) = 0. 

While this canonical prescription is completely standard in quantum field theory, it was established 
in field theories where the Hamiltonian and the fields were all Hermitian. Thus in the non-Hermitian 
but PT-symmetric case we can no longer identify the propagator this way. Rather, in the PT case one 
must distinguish between left and right vacua. Specifically, if we expand the scalar field in creation and 
annihilation operators according to <f>(x) = XX a (^) e * fc x + a(k)e~ lk ' x ) where the scalar field is not now 
Hermitian and a(k) is not the Hermitian conjugate of a(k), the left and right vacua will be defined as the 
states that are annihilated according to a(k)\Q^) = 0, (Qi,\d(k) = 0. In this case (Ol| will not be the 
Hermitian conjugate of |^r), and the propagator will be given not by the canonical (Q\T(<f)(x)(j)(x'))\Q) 
but by (Qi,\T((f)(x)(f)(x'))\Qii) instead. In the non-Hermitian but PT case (fij^ an d (^r| while distinct are 
related according to = (JIrJV, where V is the operator that relates a non-Hermitian Hamiltonian H 
to its adjoint according to VHV^ 1 = H* (see e.g. |12] and references therein). In this case the closure 
relation and normalization conditions are of the form 



Y^\ni)(ni\V = I, (n t \V\n,) = 5 tv (12) 

i 

and the propagator is given by (f^ \T((f>(x)(f>(x')) = (Q^\VT((f)(x)(j)(x'))\fl^) . With there being no 
states with (ni\V\rij) = —Sij, there are no states of negative norm, with it instead being the presence of 
the V operator that generates the relative minus sign in the propagator. Specifically, the V operator is 
related to the C operator of PT theory an operator that commutes with the Hamiltonian and obeys 
C 2 = I, with the relative minus sign in the propagator being due to the fact that the eigenvalues of the C 
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operator are +1 and —1. The relative minus sign is thus associated with the operators and not with any 
negative norm states. Thus what had led to the belief that the Ds(k 2 ) propagator had to be associated with 
negative norm states was the assumption that one had to identify the propagator as (Q\T(<f)(x)<j)(x'))\Q) 
and had to use the Dirac norm. While such a Dirac norm is appropriate in the Hermitian case, as we see, 
the non-Hermitian situation is much richer, with the bra state in a scalar product not in general needing 
to be the Hermitian conjugate of the ket. 

As can be seen from the structure of the scalar field theory action /g and propagator D$(k 2 ), not only, in 
analog to the Pais-Uhlenbeck Dp\j(u 2 ), do we encounter both a relative minus sign (when M 2 is non-zero) 
and a singularity (when we set M 2 = 0), in addition we see that in the M 2 — > limit the second-order 
term in the Jg action drops out, with the action becoming pure fourth order. Since the identification 

<t>(x,t) ~ z(t)e ihs , cut = (P + M 2 ) 1/2 , co 2 = \k\ (13) 

(and thus uj\ — uj\ = M 2 ) reduces the scalar field theory to the Pais-Uhlenbeck theory, we see that study of 
the Pais-Uhlenbeck theory enables us to study unitarity issues both in second- plus fourth-order derivative 
quantum field theories and in pure fourth-order ones. We thus turn now to a detailed study of the Pais- 
Uhlenbeck theory, and after establishing that fourth-order theories are unitary, we make an application of 
the pure fourth-order conformal gravity theory to galactic rotation curves. 



2 The Pais-Uhlenbeck Oscillator Model 

As introduced, the Pais-Uhlenbeck action Jpu contains three degrees of freedom, z, i, and z. This is too 
many degrees of freedom for one oscillator but not enough for two. Consequently the theory is constrained. 
To handle the constraints we replace z by a new variable x, and using the method of Dirac constraints 
obtain [TJ [2] a Hamiltonian of the form 

2 

ZJ Px , ,7/2, 2\ 2 7222 /i A \ 

Hpxj = — + p z x + - [uj 1 + u 2 ) x - -u x u 2 z , (14) 
with two canonical pairs that obey 

[x,p x ] = i, [z,p z ] = i. (15) 
If we proceed canonically we can introduce two sets of creation and annihilation operators 

z = ai e~ iwit + a\e iwit + a 2 e- iU2t + a\e iU2 \ 
Pz = i-ftu^ia^- 1 ^ - a{e iuJlt ) + i^lu^e-^ 1 - a\e iW2t ), 
x = -iuj x {axe- iu}xt - a\e iullt ) - iu 2 (a 2 e- iuJ2t - a\e iW2t ), 

p x = ^^(axe-^' + ale^O-T^^e-^' + a^^ 2 *) (16) 
to find that the Hamiltonian and commutation relations take the form 

Hp V = 27(w 2 - ^(ufalax - u^a^) + (o>i + u; 2 )/2, (17) 

[ai,a\} = - -\ 2-, [a2,4] = _ o 7~2 2V ( 18 ) 



27CU1 {UJ 2 — Cl> 2 ) ' ~" 27^2 (w 2 



IjJ- 



2, 



While we have thus succeeded in diagonalizing the unequal-frequency Hamiltonian via second quanti- 
zation, we see that the limit of equal frequencies is singular. We will thus need to treat the equal-frequency 
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case separately, and shall thus discuss the unequal- frequency case first. For the unequal-frequency case we 
note the relative minus sign between the two commutators. Thus in order to avoid negative norm states 
we must unconventionally define the vacuum as the one that a\ and a\ annihilate (we take oj± > U2 for 
definitiveness), with the energy and norm of the vacuum then being given by 



<n|fi) = a\\Q) = 0, (Q\4a 2 \Q) = > 0, # PU |fi) = ~{u x - u 2 )\Q). (19) 



In such a vacuum states of the form (<22) n |fi) have negative energy, with the energy spectrum being 
unbounded from below, just as is to be anticipated given the presence of the — {^/2)u\u 2 z 2 term in Hpy. 

To avoid such an outcome we could instead work in a Hilbert space in which a\ and a 2 annihilate the 
vacuum. However, then we obtain 



a x \n) = a 2 \Q) = 0, (Q\a 2 al\Q) = f ^ < 0, H PV \Q) = h Ul + u 2 )\Q). (20) 



While we no longer have any states with negative energy, we instead now have states with negative norm. 

As we see, both of the above realizations of the commutation algebra are undesirable and the problem 
even appears to be insurmountable. However, quantum mechanics is a global theory, and we must thus 
take global considerations into account. To this end we make a Schrodinger wave mechanics representation 
of the commutation algebra of the form 

■ 9 .0 

to obtain a stationary state with energy {u)\ +u 2 )/2, viz. the ground state in the realization of the theory 
where energies are bounded from below, with its wave function being given by ip (z,x)e~ i{uJl+UJ2)t/2 where 



i/jo(z, x) = exp 



^(tUi + u 2 )uiu 2 z 2 + i^UiUJ2ZX - + U2)x 2 



(22) 



As we see, while this wave function has good asymptotic behavior at large x, it diverges at large z, and is 
thus not normalizable. Since we can only represent the momentum operator as — id z when it acts on wave 
functions that are well-behaved at large z, we see that this representation for p z is invalid. In an integration 
by parts of / dzip* (z)(—id z )ij)(z) we would not be able to throw away surface terms, and thus have to con- 
clude that in the bounded from below energy sector the operator p z cannot be Hermitian. Similarly, in this 
sector the vacuum normalization evaluates to (fl\Q) = (Q\ J dzdx\z,x)(z, x\Q) = J dzdxip^z, x)ipo(z , x) , to 
thus be infinite. Thus the global information that is missing from ( 120]) is that we cannot normalize 
to a finite value, and that accordingly no closure relation of the form given in ( II ip can hold. 

While ipo{ z , x ) is n °t normalizable when z is real, we note that it would be nicely bounded if z were pure 
imaginary. Hence the theory will be well-behaved if take p z not to be Hermitian but to be anti- Hermitian 
instead. With the z operator then being anti-Hermitian too, the troublesome —{^j c X) i J\ijj\z 2 term now is 
bounded from below. Since we can continue to maintain Hermiticity for x and p x , we see that because 
of the p z x cross term Hpu is not Hermitian. However, its eigenvalues are real, so Hpu emerges as a PT- 
symmetric Hamiltonian. Since we had made the creation and annihilation operator expansion of ( fl6|) on 
the assumption that z and p z were Hermitian, we see that ( TT6|) is not a valid expansion in the realization 
of the theory in which the energy is bounded from below. Since use of ( I16p led us to states with negative 
norm, we must seek an alternate formulation of the theory, one that uses some entirely different norm. 
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To find such a norm it is convenient to work not with anti-Hermitian z and p z but to instead introduce 
Hermitian operators y and q according to 



y = e np z z/2 ze -n Pz z/2 = _^ g = ^ Pz z 12^-^/2 = ^ (33) 

These operators also form a canonical pair with [y, q] = i, with the transformation in (1231) being one of a 
broad class of transformations that leave quantum commutators (and equally classical Poisson brackets) 
unchanged. Specifically, a commutator such as [z, —id z ] is left unchanged if z is transformed to e %e z since 
d z is then transformed to e~ l0 d z . In and of itself such a transformation would be without content unless 
in varying 9 we can transit between regions in the complex z plane with differing asymptotic behavior for 
the wave functions. The complex z plane thus breaks up into regions known as Stokes wedges [H] with 
differing asymptotic behaviors, and one is able to represent the [z,p z ] = i commutator in the differential 
form [z, —id z ] only when it acts on states ip(e t0 z) in a domain in the complex z plane where ip(e l9 z) is 
asymptotically bounded. 

In terms of the operators y and q and an x and a p x = p that obey [x,p] = i, the Hamiltonian takes 
the non-Hermitian form 

2 

H = |- - iqx + | (cof + coj) x 2 + ^cofcoly 2 , (24) 

and its PT symmetry is made manifest via the assignments: p and x are odd under P, q and y are P 
even, p and y are T odd, q and x are T even. Finally, using the techniques of PT theory we introduce a 
Hermitian operator Q defined as [6] 

Q = a[pq + -f 2 u 2 u%xy}, a = — - — log f UJl + u A (25) 
to find that the unequal-frequency Hamiltonian H given as 

H = e-^ 2 H PV e^ 2 = f + + + (26) 

27 27cuf 2 2 

is Hermitian and diagonal, with all of its energy eigenvalues being real and positive. The unequal-frequency 
Hpu while not itself Hermitian can be brought to a Hermitian form by a similarity transform. 

Since Hermitian theories possess no states with negative norm, it follows that the unequal-frequency 
Hpu theory must possess none either. To see what form the -f/pu theory norms do take, we note that since Q 
is Hermitian rather than anti-Hermitian, the unequal frequency Hpu has been diagonalized by a similarity 
transformation that is non-unitary. Since similarity transformations do not preserve orthonormality of 
basis vectors, the eigenstates of H and H are not unitarily equivalent. Rather, if we define eigenstates of 
Hpu and H according to Hpu\n) = E n \n) and H\h) = E n \h), we find that these states and their conjugates 
are related according to 

H\n) = E n \n), H PlJ \n) = E n \n), |n) = e Q/2 \n) 

(n\H = E n (n\ } (n| = (n\e Q/2 , (n\e~ Q H Pl j = (n\e~ Q E n (27) 

with the Hermitian conjugate (n\ of the right-eigenstate |n) of Hpu not being a left-eigenstate of Hp\j. 
In terms of the operator V we introduced earlier we can identify V = e~®, and can check directly that 
VHpuV^ 1 = Hp v just as required. (In passing we note that in this case the C operator of PT theory can 
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be written in terms of Q and the parity operator P as C = e®P [5J, with it obeying C 2 = I since Q obeys 
PQ = —QP.) Finally, since H is unitary, we obtain normalization and closure relations of the form 

(n\rh) = 5 m<n , E\n){n\=I, H = E\n)E n (n\, 

(n\e- Q \m) = 8 m , n , Z\n}(n\ e - Q = /, H Pl] = E\n)E n (n\e- Q (28) 

with there being no states of negative norm, and with the similarity transformation thus taking us between 
a skew basis and an orthonormal one. In the PT realization of the unequal-frequency theory then there are 
no states of negative energy and no states of negative norm, and the theory is fully physically acceptable. 



3 The Singular Equal-Frequency Limit 

Inspection of (|2"5|) shows that in equal-frequency limit the operator Q becomes singular, with the trans- 
formation to H no longer being valid. To see what happens in the limit we consider the two one-particle 
states of the realization of the unequal frequency theory in which the energy is bounded from below, viz. 



V>i(z,y,i) = (x + u 2 y)ip (x,y,t)e lu)1 , E 1 = E + ui 



ifa(x, y, t) = {x + u iy )Mx, V, t)e-^\ E 2 = E + u 2 , (29) 
where if) (x,y,t) is the ground state wave function 



fk(x,y,t) = exp 



--{ux + u^x -ju x u 2 yx - -(wi + u 2 )u 1 u 2 y - iE t 



(30) 



with energy Eq = (ui + u> 2 )/2. As we see, in the limit not only do the two one-particle energies become 
equal, the two wave functions become identical with the two one-particle states having collapsed into 
a single state. To confirm that the Hamiltonian has thus become a non-diagonalizable, Jordan-block 
Hamiltonian, we note that in the unequal-frequency one-particle sector the Hamiltonian takes the form [7] 

#lPl e ) = 7^- .2 A , ,2 Jl ), (31) 



where we have set uj\ = u + e, u 2 = u — e. For this Hamiltonian the eigenstates are given by 

|2. + £ > = ( 2 7 f ), |2.- f >.( 2 ^ £ ) (32) 

corresponding to eigenenergies 2u + e and 2u — e. 
If we introduce the operator 

CM - 1 f2u + e -(4u 2 -e 2 )e\ 

b[€) ~ 26^/2(20; + e)V2 v e (2u + e)e 2 ) {66) 

we find that i?ip(e) is diagonalized by this S(e) according to 
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Then, in the e — > limit, we see that S(e) becomes undefined and H 1P (e) becomes the Jordan-block 
Hamiltonian 



H 1P (e = 0) = 2co 



(35) 



with two degenerate energies but only one eigenstate. (As noted in j7], by setting y(e = 0) = e~ lujt [— i(a — 
b) — 2but] + e luJt [— i(a — b)+ 2bu)t] one can introduce creation and annihilation operators at e = that obey 
[a, a] = 0, [b, b] = 0, [a, b] = [b, a] = I/87W 3 , to then obtain Hp V (e = 0) = 8'yu 4 (2bb + ab + ba) + u, with 
this form for Hpu(e = 0) updating the form given in |TJ since a and b are not the Hermitian conjugates of 
a and b. With this form for H P u(e = 0) the one-particle sector i?ip(e = 0) then follows.) 

To ascertain what happened to the eigenstate of H\p{e) that went missing in the limit, we take the 
limit of the ground state and of two particular linear combinations of the two one-particle eigenfunctions 
of the unequal-frequency theory, to obtain 



lim?po(x,y,t) = exp 

e— y\j 



—•yuj 3 y 2 — ^u^yx — •yux 2 



iut 



(36) 



^ia(x,y,t) 



hm = {x + uym\n,y,t)e 

e— 5-0 Z 



lim 



ifj 2 {x,y,t) - V>!(a;,y,t) 
2e 



[(x + uy)it + y] ip (x,y,t)e 



-iwt 



(37) 



As we see, the non-singular combination ipi(x,y,t) remains an eigenstate in the limit, but the singular 
combination ipi a (x,y,t) does not. However, while i[)i a (x,y,t) is not an energy eigenstate, it is a non- 
stationary solution to the time-dependent Schrodinger equation 



Q „ fid 2 d 7 \ - 

z-VOr, y, t) = (-^^2 ~ X dy- + 7wV + 2 ^ ) V " t) ' 



(38) 



Despite the fact that the non-stationary states are not energy eigenstates, the set of stationary plus non- 
stationary states combined is complete since together they provide just the right number of independent 
polynomial functions of x and y that are needed to form an arbitrary localized wave packet in the (x, y) 
plane [7]. Moreover, as shown in [7J, the norm of any such packet is preserved in time, with time evolution 
being unitary. As we thus see, in the event that a Hamiltonian is not diagonalizable, unitarity is still 
expressed as a completeness statement. However, now it is the set of solutions to the time-dependent plus 
time-independent Schrodinger equations combined that needs to be complete rather than solutions to the 
time-independent one alone. 



4 Conformal Gravity 

The conformal Weyl tensor defined as 

is a quite remarkable geometric object. Under a local conformal transformation of the form 

g„u(x) -> e 2a ^ x) gfMU (x) 



(39) 



(40) 
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with arbitrary spacetime-dependent a(x) the Weyl tensor transforms as C\^ UK — > e 2a ^C\^ UK . Thus even 
though the Riemann tensor R\^ UK and its R^ K = g Xu R\^ U K and R a a = g^ K R^ K Ricci contractions all acquire 
derivatives of a(x) under the conformal transformation (each being a second-derivative function of the 
metric), all of these derivative terms cancel identically in the Weyl tensor. As such, the Weyl tensor 
behaves under a conformal transformation the same way that the Maxwell tensor F^ v behaves under a 
gauge transformation. It is thus very suggestive to endow gravity with such a conformal structure, so that 
as well as being general coordinate invariant, gravity will in addition possess local conformal invariance as 
well. With such additional symmetry the gravitational action is uniquely prescribed to be of the form 

iw = -a g J d A x{-gfl 2 C HlVK C x ^ = -2a g J d^-g) 1 ' 2 

where the gravitational coupling constant a g is dimensionless, with no other terms being permitted. (To ob- 
tain a dimensionless a g one needs four derivatives to compensate for the / d A x measure, with the Cx fJiUK C x>J,UK 
term precisely providing the needed four derivatives). Associated with the Jw action are fourth-order 
derivative gravitational equations of motion of the form (see e.g. |4]) 

Aa g [2C^ XuK . x . K - C^Rxn] = 4a g W IJV = T" u , (42) 

where W^ v is covariantly conserved and traceless and can be written in the form 

- jj<n# a J ; % + l(R a a )' KV + \R a a R" u - \g» v {R\) 2 - (43) 

As constructed, conformal gravity is thus a power- counting renormalizable theory of gravity since a g 
is dimensionless. Moreover, conformal gravity is able to control the cosmological constant since conformal 
invariance does not permit the presence of any I = J d 4 x(— g) 1//2 A term in the fundamental action. In 
addition, we note that since the cosmologically relevant Robertson- Walker (RW) geometry is a conformally 
flat geometry (viz. one in which C\^ VK = 0), in conformal cosmology (1421) reduces to 

T>» = 0, (44) 

so unlike the (non-conformal) double-well Higgs potential, conformal gravity knows exactly where the zero 
of energy is, a feature of the theory that is central to solving the cosmological constant problem (T3J [5] . 

As well as exclude any cosmological constant term from the fundamental action, the same conformal 
symmetry also excludes the presence of any Einstein-Hilbert term /eh = —(c 3 /IQitG) J d 4 x(—g) 1 ^ 2 R a a in 
the fundamental action. Despite this, the theory still admits of the Schwarzschild solution jHJ [15] since 
Rp,v = is an exact solution to fj42l when T^ v = 0. The conformal theory thus recovers both Newton's 
Law of Gravity and its familiar general- relativistic corrections (gravitational bending of light, gravitational 
redshift, and modifications to planetary orbits). As noted in [4], Einstein gravity is thus only sufficient to 
give the standard Schwarzschild metric phenomenology but not necessary. 

Now while the vanishing of R^ u implies the vanishing of W^, the converse is not true since the vanishing 
of W^ u can be achieved without R^ v vanishing. The vacuum vanishing of W^ v thus has other solutions 
[141 [15], solutions that can thus enable us to distinguish between Einstein gravity and conformal gravity. 
As we shall see below, these other solutions are important on galactic distance scales where they enable 
the conformal theory to account for galactic rotation curve systematics without the need for dark matter. 



[IV 



3^ 



(41) 
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To identify the quantum structure of the conformal theory, it is convenient to linearize the theory 
around flat spacetime according to g^ v = rj^ + h^ v where r]^ u is the Minkowski metric and h^ u is the 
fluctuation. Through second order in h^ u it is found [16] that (jHJ) and ( jUJ) only depend on the traceless 
combination K^ v = h^ v — {l/4)ri ilv T)' x °h( X p- The equations of motion are found to simplify considerably in 
the transverse gauge where d^K^ = 0, with the first order modification to W^ u being given by 

W^{1) = l -{d a d a ) 2 K» v (45) 

(with both sides of (I45p being traceless), and with the second-order fluctuation term in the conformal 
action Jw being given by 

J w (2) = - y / d A xd a d a K, v d^K^. (46) 

As we see, both fj4"5]) and f l4"6]) are diagonal in the (/i, v) indices, and so each component of K^ u can be 
treated independently. Comparing with ([7]) and (IE]), we see that each independent component of K^ u obeys 
the same dynamics as a fourth-order scalar field with M 2 = 0. The unitarity analysis of the scalar field 
theory this caries over identically to the conformal theory [5]. Thus despite the fourth-order nature of its 
equations of motion, the unitarity of microscopic conformal gravity is secured. Since the macroscopic limit 
of the conformal theory is obtained by taking matrix elements of its quantum field operators in states with 
an indefinite number of gravitational quanta, the consistency of conformal gravity as a macroscopic theory 
is thus secured as well, and so we turn now to a study of some of its macroscopic implications. 



5 Fitting Galactic Rotation Curves 

To fit galactic rotation curves one needs to determine the geometry for a static, spherically symmetric 
system in the conformal case. To this end it was noted in [T4] that though use of the underlying conformal 
symmetry the most general metric in the static, spherically symmetric case can be brought to the form 

ds 2 = B(r)dt 2 - ^— - r 2 d0 2 - r 2 sin 2 6d<j) 2 , (47) 
B(r) 

to thus be described in terms of just the one function B(r) where r is the radial coordinate. Moreover, in 
[15] it was shown that in metric given in (T4T1) the gravitational equation of motion given in ( 1421) takes the 
form of a fourth-order Poisson equation 

' (W\-W r r ) = B"" + = V 4 5 = -— — (T° -T;) = /(r). (48) 



B{ry u r> r 4ot g B(r, 

As derived, f T4Hl) is exact and without approximation, showing that despite its complexity, in the presence 
of sufficient symmetry (T42]) can be quite tractable. Moreover, for a general source function f(r) (148]) can 
be integrated exactly to yield 

B(r) = — / dr'r' 2 f(r') / dr'r' A f{r') - - / dr'r' 3 f(r') / dr'r 1 f \r') + B (r) , (49) 



2 jo 6r Jo 2 Jr 6 

dr'r' 2 f{r') + — df V 4 f \r') - - dr'r'f(r') + B'(r), (50) 



dr 2 Jo 6r 2 jo 3 
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where B(r) satisfies V 4 B = 0. 

In contrast to (jUJ) and (15 Op . we recall that for the second-order Poisson equation V 2 0(r) = g(r), the 
solution is of the form 

0(r) = — / dr'r' 2 g(r') - / dr'r'gir'), (51) 

r JO Jr 



d(p{r) 



\ f dr'r' 2 g(r'). (52) 
r A Jo 



dr r 2 Jo 

Comparing fl50l) and fl52|) . and noting the presence of an integral that continues all the way to infinity 
in fl50l) . we see that the novel feature of the fourth-order theory is that the gravitational force at radius 
r receives contributions not just from local sources in the < r' < r region but from global sources in 
the r < r' < oo region as well. Thus unlike Newtonian gravity where local gravity is produced solely by 
local sources, in the fourth-order case global sources can contribute as well. Newtonian gravity is thus 
intrinsically local, while conformal gravity is intrinsically global. Since in Newtonian gravity one can only 
explain local gravitational observations though the use of local sources, a failure to account for observations 
using luminous sources alone leads to the need to postulate the presence of dark matter within any given 
system of interest. However, in the conformal theory there are additional contributions to the gravitational 
field within systems due to matter that is outside of them. As shown in [T7J dHl [TjJJ [2D], these global effects 
can remove the need for galactic dark matter, with dark matter potentially being nothing more than an 
attempt to describe global physics effects in purely local terms. To determine the effect of conformal 
gravity on galactic motions we thus need to determine the local contribution due to the material within 
each given galaxy and the global contribution due to the material in the rest of the Universe. 

As regards first the local contribution, we note that if in ( j4"9|) we consider a purely localized source of 
radius tq, then on setting B(r) = 1, in the region r > tq we obtain 

28 

B(r > r ) = 1 - — + 7T, (53) 
r 

where 

1 rro 



7 = "2 



/ dr'r' 2 f{r'), 2(5 = - / dr'r' A f{r'). (54) 
Jo 6 Jo 

In (|53p we obtain a Newtonian potential term even though there is no reference to any second-order Poisson 

equation. Thus, as noted earlier, standard second-order gravity is only sufficient to give the Newtonian 

potential but not necessary. In addition we obtain a linear potential term. This term modifies standard 

gravity at large distances but not at small ones, i.e. it modifies gravity in precisely the region where the 

standard theory has to resort to dark matter. 

In conformal gravity a local source such as a star produces a potential of the form 

= ( B( r) - iy = _r* -av 

K ' 2 r 2 K ' 

for each solar unit of material. For spiral galaxies the luminous material is typically distributed in a disk 
with surface brightness = Soe"^^ at radial coordinate R. Integrating the potential of ( |55l) over 

such a distribution yields a contribution to circular velocities in the plane of the disk of the form [1] 



^LOC(^) 



N*(3*c 2 R 2 
2R 3 o 



2Rq J V 2Rq J V 2Rq J V 2R\ 



W. h (JL) Kl (JL , , m 



2Rci \2RqJ \2Rq 
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where N* = 2tt(M/L)'E Rq/Mq is the total number of stars in a galaxy, M/L is the galactic mass to light 
ratio, and M = N*M & is the visible galactic mass. 

As regards the global contribution, we note that it consists of two components, a contribution from the 
homogeneous background Hubble flow, and a contribution due to inhomogeneities in it such as clusters of 
galaxies. Since W^ v vanishes identically in an RW geometry, the Hubble flow only contributes to the B{r) 
term in (|49[) . To determine the effect of the Hubble flow on motions within galaxies we need to transform 
from comoving Hubble flow coordinates to coordinates at rest in a galaxy. With the transformation 

P = oTi I \T7TT~n~, > r = / dtR(t) (57) 



2(1 + 7o r) 1 /2 + 2 + 7o r' 



effecting 



,2 dr 2 „ 1 /l + 7o p/4\ 



2 r 



(1 + 7 r)or - r - r z dQ 2 



dr< - fl2 2 ( aL 12 (dp 2 + p 2 dn 2 
[1 - 7oP /W V 



(5? 



(l + 7o r) R 2 (r) \l- l0 p/4j 

we see that in the rest frame a topologically open RW cosmology with negative 3-curvature K (viz. the one 
preferred in conformal cosmology [I]) acts just like a universal linear potential with strength 7 = 2(— K) 1 ^ 2 . 

The inhomogeneous contribution is generated by the fourth integral in ([4*9]) . And in the region r < r c \ us 
where r c i us is a typical cluster scale at which the global f(r') begins, this integral acts as a quadratic 
potential — nr 2 with universal strength k = (1/6) dr'r' f(r'). Finally, adding the contributions of the 
global linear and quadratic potentials to the local (1561) yields for the total velocity and its asymptotic limit 

v 2 OT (R) = v 2 LOC (R) + ^ - kc 2 R 2 , (59) 

2 N*p*c 2 N*Yc 2 R loc 2 R 2 2 
v TOT (R) -> — ^— + + — kc R . (60) 

Despite the fact that the galaxy- dependent N* is the only free parameter in 059|) . using the values 

f3* = 1.48 x 10 5 cm, 7* = 5.42 x RT^cm -1 , 70 = 3.06 x lO^cm -1 , « = 9.54 x 10~ 54 cm- 2 (61) 

Mannheim and O'Brien [TDl |2"U] were able to fit the rotation curves of a large and varied sample of 
138 spiral galaxies using only the visible galactic N* and no dark matter at all, viz. just one parameter 
per galaxy fits. In Fig. (1) we show fits to 21 of the largest galaxies in the 138 galaxy sample. Since the 
ft-dependent term appears with a negative sign in ( 159]) and (160"]) . there would have to be an iV*-dependent 
bound on the size of a galaxy with iV* stars since i4ot(-^) cannot be negative. To illustrate this effect we 
have included in Fig. (1) Mannheim and O'Brien's previously unpublished extended distance expectations 
for the three galaxies in the sample with the largest N* values. The novel prediction of an explicit end to 
galaxies provides a falsifiable signal for the conformal gravity theory. 

That we are able to fit the data at all with a formula as constrained and universal as ( 159]) is very 
encouraging for the conformal theory, and especially so since, in addition to each visible galactic mass, dark 
matter fits to the same 138 galaxy sample require an additional 276 free parameters (two free parameters 
for each galactic dark matter halo). The reason why the fits do succeed is because across the entire 
galaxy sample the value of v 2 :OT (R)/c 2 R at the last detected data point of each galaxy is numerically 
found [HI [19j [20] to not only be close to universal but to be close in magnitude to the universal 70. Such 
universality is a property of the data themselves, and has yet to be accounted for by dark matter theory. 

To conclude, we note that the fits not only provide support for the viability of conformal gravity at the 
macroscopic level per se, but in addition they provide evidence for the relevance to physics of an underlying 
microscopic theory that is based on a non-diagonalizable but PT-symmetric Hamiltonian. 
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